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ABSTRACT 



We present results of high resolution (~ 55000) spectral observations of 830 
photometrically pre-selected candidate red giants in the magnitude range of V 
= 9-12. We develop a pipeline for automated determination of the stellar atmo- 
spheric parameters from these spectra and estimate T e ff, log g, [Fe/H], micro- 
turbulence velocity, and projected rotational velocities, vsini, for the stars. The 
analysis confirms that the candidate selection procedure yielded red giants with 
very high success rate. We show that most of these stars are G and K giants with 
slightly subsolar metallicity ([Fe/H] ~ -0.3 dex). An analysis of Mg abundances 
in the sample results in consistency of the [Mg/Fe] vs [Fe/H] trend with published 
results. 

Subject headings: stars: abundances, fundamental parameters, atmospheres, late- 
type, rotation 



1. Introduction 



The photometric Grid Giant Star Survey (GGSS; [Patterson et al.ll200ll ; iMajewski et al. 



20001 ) was designed to search for G and K giant stars in selected relatively large areas (~ 
1 sq. degree), more or less evenly spaced all around the sky. A number of GGSS giants 
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were selected as potential candidates to be astrometric reference grid stars for the NASA 
Space Interferometry Mission (SIM). The primary motivation for the spectroscopic follow-up 
survey of GGSS candidate giants was to monitor the line-of-sight velocity at a high level of 
precision, using an iodine cell to identify radial velocity - stable stars in support of SIM. 
During the radial velocity monitoring campaign of the Northern hemisphere sample at the 
McDonald Observatory (University of Texas), we obtained some 3000 spectra for more than 
1500 candidates with magnitudes V ~ 9 - 12. Although an iodine cell was used to imprint 
an I 2 spectrum onto the stellar spectrum in order to monitor radial velocities, a significant 
number of spectra were obtained without the cell; these have moderate signal-to- noise (S/N), 
and are suitable for an abundance analysis. The spectroscopic analysis of this large number 
of spectra is the object of this paper which discusses an automated analysis for deriving 
stellar parameters and metallicities. 

We have previously analyzed, using an automated tec hnique, the ste l lar pa rameters for 
an earlier sample of SIM reference grid star candidates (IBizyaev et al.l 120061 ) . The stars 
considered in the present paper have somewhat brighter apparent magnitudes, because the 
selection criteria were changed to accommodate revised SIM technical requirements, and 
in gen eral have higher S/N (typically S/N ~ 30 - 50). Also, in contrast to our earlier 
study (IBizyaev et al.l 120061 ). we employ here a different and im proved technique to derive 
the parameters using an approach similar to that developed by (IKatz et al.lll998l ). which is 
found to be better applicable to spectra with S/N of moderate quality (of the order of tens). 
In this paper, we focus on the characterization of the GGSS sample with the construction 
of a pipeline for the derivation of the target stars atmospheric parameters and metallicities 
(iron abundances). Such results will add to the large database of metallicities which have 
been produced oy er recent years by large surv eys using low-resolution spectra (e.g. SDSS) 
( iLee et al.l l2008al Jbl; lAllende Prieto et al.ll2008l ). The high-resolution spectra in this study, 
in contrast with surveys which analyze lower resolution spectra, will allow future studies to 
determine the abundances of several chemical elements with accessible spectral lines in the 
wavelength range covered by these spectra. 



2. Observations 



For the preparatory RV campaign for SIM, we conducted spectroscopic observations of 
a large sample (1500+ objects) of photometrically pre-s elected candidate red giants. The 



whol e-sky sample is a mix of the brightest GGSS stars (iFrink et a 
2001 ). as well as Tycho stars ( Catanzarite et al. 2004 ; Ciardi 12004 ; 



2001 



Patterson et al. 



Thompson et al.ll2005l ). 



Note that this joint sample has only a very small overlap with that studied by Bizyaev et 
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al. (2006) because of the revised SIM technical requirement to use brighter reference stars. 

Observations were conducted using the 2.1m telescope at the McDonald Observatory 
and the Sandiford Cassegrain echelle spectrograph, which provides a nominal resolution R 
= 55000 ( McCarthy et al. 19 93). In contrast to the preceding radial velocity monitoring 



campaign (IBizvaev et a 



nique (e.g. 



Butler et al. 



. 2006), the observations discussed here employed the iodine cell tech 



19961 ) . The gas cell (provided by W. Cochran, Cochran fc Hatzes 



1994 ) with evaporated iodine was heated to +50°C and placed in the beam before the slit 
of the spectrograph. For most of the stars, observations were made with and without the 
iodine cell. The latter frames were used as templates for the radial velocity measurements. 
The radial velocity monitoring results will be the subject of a future paper. Here we focus 
on the echelle spectra obtained without the iodine cell and which will be used for abundance 
analyses. We observed 841 spectra for 830 program stars. In addition, we obtained spectra 
of 10 very bright stars with previously identified atmospheric parameters for comparison 
purposes. A day-light solar spectrum is also included in the analysis. 

Th-Ar comparison spectra were taken right before each program star. The adopted 
instrumental setup enabled us to cover the spectral range from 5000-5900 A . There were 
24 spectral orders across the chip, each order being 1200 pixels long. The first and last (as 
a rule incomplete) orders were rejected, as well as three orders in the middle of the CCD 
chip that are contaminated by an internal reflection. The mean S/N for all spectra was 42, 
with an r.m.s. of 12, and the typical exposure time was 1000 sec. Observations of the bright 
comparison stars resulted in spectra with S/N of 150 or more. 

The observed spectra were reduced from two-dimensional to "echelle" format in a stan- 
dard way with the IRAF0 software package. Corrections for bias, flat field, and scattered 
light were applied, and cosmic ray hits cleaned out by the IRAF's tasks from crutil package. 
The internal accuracy of the wavelength calibration via the Th-Ar reference spectra was on 
average 0.001 A. 

Intermediate resampling introduces systematic errors in the wavelength calibration and 
degrades the RV accuracy. We preserve the initial sampling and do not concatenate the 
echelle orders into a single 1-D spectrum. Rather than that, we keep each echelle order 
separate, with its individual sampling. 



1 IRAF is distributed by NOAO which is operated by AURA, Inc. under contract with the NSF. 
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3. The TAGRA Pipeline 

We designed an automated pipeline to derive the effective Temperatures, iron Abun- 
dances, and surface GRAvities (TAGRA) for the sample stars via comparisons of observed 
spectra with synthetic ones. Near-IR photometric data were also used in the pipeline (see 
below) . 

Observed spectra were normalized to an empirically defined continuum le vel. We used 



the ul tra-high resolution spectra of both the Sun and Arcturus from the atlas by lHinkle et al. 



(120001 ). degraded to our resolution, in order to identify pieces of spectra relatively free of 
spectral lines and which could be used as continuum points. For each spectral order we used 
as many continuum points as possible; Table 1 lists 235 such continuum points which were 
included in the pipeline. The continuum points in our observed spectra were fit with low- 
order polynomial functions. Figure [1] shows an example of the continuum fitting within one 
spectral order of the degraded Arcturus spectrum, whose spectrum is quite representative of 
the stars analyzed here. 



3.1. Effective Temperatures 



The introduction of one more degree of freedom (T e fj in this case) into the modeling 
and fitting using synthetic spectra leads to larger uncertaintie s in the derived stell ar pa- 
rameters. Since all program stars were observed by 2MASS ( ISkrutskie et al.ll2006l ). this 
provides additional photometric information which can be used to const rain effective tem- 
peratures for the pipeline. Similar to the analysis by iBizyaev et al.l (120061 ) . effe ctive temper 
atures were determined from the 2MASS colors (J-K). We use a calibration by lAlonso et al. 
(199 9|) and transformation equations from 2MASS to Alonso's TCS color system given by 
Ramirez fe Melendezl (120051) . The galactic foreground extin ction is taken from the GGSS and 
Tycho sample papers ( jPatterson et al.l l2001t ICiardil I2004J ). We note that this calibrations 
of (J-K) for giants does not have a dependence on the metallicity. Since the program stars 
are relatively bright, this calibration provides good accuracy for the effective temperature, 
typically better than 125 K (Alonso et al. 1999). 



3.2. Iron Lines, Surface Gravities and Microturbulent Velocities 



With effective temperatures defined by the 2MASS colors, iron lines (Fe I and Fe II) 
were used to determine the remaining stellar parameters: surfa ce gravity and microtu rbu- 
lent velocities. Potential Fe lines were taken from the studies by iFulbright et al.l (120061 ) and 
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Bizyaev et al.l ( 120061 ). The spectral regions containing these Fe I and Fe II lines were in- 
vestigating for their suitability as Fe abundance indicators by using the complete linelists 
(containing all lines from all species from Kurucz fe Bell 1995) and generating synthetic 
spectra with the M OOG spectrum synthesis code (jSnedenl Il973l ) and ODFNEW Kurucz's 
atmospheric models (jCastelli fe Kuruca 120031 ). A small spectral window around each Fe line 
was studied as to how the line responded to variations in the surface gravity, iron abundance, 
and microturbulence (log g, [Fe/H], and £). Spectra in each of the selected regions were 
inspected and any heavily blended Fe lines were excluded. Also very weak lines that were 
not well-defined in our medium quality spectra were excluded from the final analysis. We 
ended up with 81 spectral regions; the list of wavelength ranges and the Fe lines in each of 
them is presented in T able 2. The excitation pot entials and logfgf ) valu es for the sample Fe 
lines are adopted from iBizyaev et al.l (120061 ) and iFulbright et al.l (120061 ) . 



Radial velocities are determined as part of the spectral analysis and initial measurements 
of the stellar radial velocities were determined from synthetic spectra calculated for the region 
around the Mg I 5180A line. The synthetic spectra were then convolved with a Gaussian 
instrumental profile, which corresponded to the typical FWHM of the observations (2.4 
pixels); the radial velocity (RV) of the star was then determined via cross correlation of 
synthetic and observed spectra. This RV was used as the initial radial velocity in further 
modeling. 

As the next step in the analysis of a given star, synthetic spectra were computed for each 
spectral window in Table 2 and were fit to the observed spectra using as free parameters RV, 
log g, [Fe/H], £, and vsini. Note that since the echelle orders have overlapping wavelengths, 
a few sample lines were observed in more than one order. To avoid using low S/N regions 
that can occur near the edges of echelle orders, we rejected fitting a line if it was closer than 
100 pixels to the edge. Variations in the width of the instrumental profile (IP hereafter) 
were also investigated as an additional model parameter. Because the spectrograph focus 
was monitored regularly during the observations with Th-Ar line profiles, we found that the 
IP was stable enough during observations, so that it could be kept fixed for each fit. We 
tried the Powell, Levenbe rg-Marquardt, as well as the downhill-simplex (also called Amoeba) 
optimization algorithms (IPress et al.lll992l ) to derive the best-fit parameters. The optimized 
values are the chi-square derived from the synthetic and observed spectra for each considered 
line region. Tests indicated that, although being the slowest of the three algorithms tested, 
the downhill- simplex algorithm yielded realistic best-fit parameters, whereas the Powell and 
Levenberg-Marquardt algorithms often failed to find the minimum. 

We considered two approaches for this analysis: 1) doing the interpolation between 
the model atmospheres to generate spectra for each set of parameters versus 2) doing the 
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Table 1. Normalized pseudo-continuum level between 5000 and 6000 A. 



X(start) 


A (end) 


Mean Cont.(Arcturus and Sun) 


r.m.s. 


4993.10 


4993.20 


0.99240 


0.00099 


5003.50 


5003.60 


0.98205 


0.01336 


5009.00 


5009.10 


0.97865 


0.00771 


5015.60 


5015.75 


0.97645 


0.00474 


5018.90 


5019.00 


0.98005 


0.01690 



Note. - The wavelength range (in A) and the normalized 
continuum level in the spectra of Arcturus and the Sun. Table 1 
is published in its entirety in the electronic edition. 



Table 2. The spectral window around the iron lines utilized in TAGRA. 



X(start) 


X(end) 


line 


A, A 


5054.55 


5054.90 


Fel 


5054.643 


5141.45 


5142.05 


Fel 


5141.739 


5236.00 


5236.30 


Fel 


5236.202 


5307.00 


5307.70 


Fel 


5307.361 


5412.65 


5412.90 


Fel 


5412.784 



Note. — Table 2 is published in its 
entirety in the electronic edition. 
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interpolation between synthetic spectra generated with a fixed grid of parameters. The 
latter approach is preferable for our computations as it is much faster and requires less 
computer resourses since we use a limited number of short pieces of spectra. In both cases 
we synthesized pieces of spectra in the selected wavelength reg ion. The syn thetic spectra 



were computed using the LTE spectrum synthesis code MOOG (jSnedenlll973l ). 



For the first method, synthe tic spectra were generat ed via an interpolation between 
the Kurucz atmospheric models JCastelli fc Kuruczl I2OO3I ) with parameters varied until a 



best fit between synthetic and observed spectra was obtained. In the second approach, 
we produced an entire grid of synthetic spectra in the selected wavelength ranges for T e ff 
=3500 - 6000 K, log g = - 5 dex, [Fe/H] = -2.5 - +1.0 dex, and £ = 0.5 - 3.5 km/s 
with steps of 250 K, 0.5 dex, 0.5 dex, and 1.0 km/s, respectively for Teff, log g, metallicity 
and microturbulent velocity. Real spectra were then interpolated between the synthetic 
spectra and these results were compared to those obtained directly from synthesis. If a 
linear interpolation was applied, the difference between the spectra generated using the two 
approaches could be as close as 1%. Using a second-order polynomial interpolation, the 
difference becomes less than 0.5%, and thus we expect to see no significant difference if 
spectra from both methods are compared with real spectra of intermediate to good quality 
(S/N up to 200). We conclude that use of the second-order polynomial interpolation yields 
accuracies sufficient for values of the S/N up to about 200. We used the second-order 
polynomial interpolation between the pre-compiled spectra in the final pipeline. 

Projected rotational velocities, vsini, were measured from fits to all spectral lines falling 
within the echelle order near A5740A (with rejection of 100 pixels from the edges of the order). 
A synthetic spectrum was generated using the best-fit model parameters for each star and this 
spectrum is broadened over a range of values of vsini), with the best fit used to define vsini. 
We note that at the spectral resolution of R=55000, this broadening which is fit as rotational 
will include underlying macroturbulence within the star as well, thus any values of vsini <8 
km-s _1 should be considered as upper limits. At the same time, the rotational velocities for 
the Sun and Arcturus derived by TAGRA are very close to their real values. Our attempts 
to fit the spectra with the rotatio nal and mac roturbulence broadening simultaneously, both 



directly and in Fourier space (see iGray Ill976l ). failed because of insufficient S/N and a lack 



of clearly isolated lines in the spectra of our stars. An alternative approa ch to the vsini 



estima tes using the cross-correlation technique for our stars is developed by ICarlberg et al. 



(120101 ) and results in the discovery of a few rapidly rotating stars in our sample. We postpone 
a more detailed discussion of the rotational and mactoturbulence broadening until the next 
paper in this series. 



The final atmospheric parameters for the sample stars as well as the estimated vsini are 
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given in Table 3. The resulting best fiting spectra to narrow spectral regions around selected 
iron lines for a typical star in our sample are shown in Figure [2j The squares indicate the 
observing spectra, the solid curves show the best-fitting synthetic models, and the dashed 
curves represent [Fe/H] offsets by ± 0.5 dex. 



3.3. Reliability of Derived Stellar Parameters 



The external accuracy of the derived stellar parameters in this study can be evaluated 
from independent determinations of stellar parameters (all using TAGRA) for a sub-sample 
of 11 stars which have been observed two times during the Survey. Table 4 summarizes the 
mean values of log g, [Fe/H], £, and vsini, and standard deviations, obtained for these stars. 
It can be seen that the typical accuracy is 0.24 dex and 0.05 dex for log g and [Fe/H], and 
0.09 and 0.17 km/s for £ and vsini. The internal precision of the derived parameters can be 
estimated by evaluating each parameter in all selected spectral regions by varying them one 
at a time. We found that the formal internal accuracy is less than external accuracy. 

In addition, high-quality spectra for a sample of bright and well studied stars were 
obtained under the same conditions and with the same instrument and setup as all other 
program stars studied here; these spectra had signal-to-noise ratios better than ~ 150. The 
stellar parameters derived here with TAGRA (and listed in Table |5]) can be compared with 
results from classical s pectroscopic analyse s previou sly published in the lit erature (listed 
in Table [5]): Arcturus (ISoubiran et al.l 120081 ). /3 Qph (Luck fc Challenerlll995l ), 7Dra, -uDra, 
vrHer, 11 Umi, HD160290, HD218029 (all from lMcWilliamlll99ol ). HD152812 farown et al 
19891 ). as well as the Sun. Figure E] shows the comparison between the published stellar 
parameters and those derived with TAGRA. The mean difference in log g, [Fe/H], and 
£ within the whole sample is -0.24 ± 0.40 dex, -0.10 ± 0.12 dex, and 0.07 ± 0.12 dex, 
respectively. Note that for the Sun and Arcturus in particular, the differences (TAGRA - 
published) for the log g and [Fe/H] is only -0.04 and -0.08 dex, respectively. These results 
are also shown in Table [5J 

As a further test, we also performed a "classical" spectroscopic analysis for one of the 
bright stars, /30ph, using the same list of lines as in the pipeline analysis. If a similar 
approach is adopted of using the effective temperature from the photometric calibration 
as a fixed value, the other derived parameters are: log g = 1.3 dex, [Fe/H] = +0.02 dex, 
and £ = 1.9 km-s -1 , i.e. in good agreement with those found with the TAGRA. Note that 
these parameters were determined under an assumption of published T e ff. If the effective 
temperature is also derived spectroscopically (from the excitation equilibrium of Fe I lines) 
we obtain: Teff=4630 K, log g = 1.9 dex, [Fe/H] = +0.15 dex, and £ = 1.9 km/s. 
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Table 3. The atmospheric parameters derived with TAGRA 



Name 


RA(J2000) 


DEC(J2000) 


m B 


A v 


T eff 


log g 


[Fe/H] 




vsini 


S/N 










mag 


mag 


K 


dex 


dex 


km/s 


km/s 




G0000+67.16328 


14:59:18.5 


67:34 


02 


10.80 


0.49 


4624 


2.32 


-0.53 


1.54 


2.46 


37 


G0024+61.2521 


00:15:06.6 


62:26 


51 


10.58 


0.45 


4489 


1.97 


-0.36 


1.41 


2.37 


38 


G0024+61.6685 


00:15:22.2 


62:18 


15 


10.76 


0.47 


4677 


3.84 


0.08 


0.50 


0.62 


28 


G0024+61.9144 


00:18:18.2 


62:13 


14 


10.42 


0.43 


4797 


2.34 


-0.11 


1.49 


1.85 


34 


G0112+61.6061 


00:47:22.8 


62:16 


03 


10.69 


0.47 


4645 


3.36 


-0.54 


0.94 


5.68 


14 


G0121-16.44 


00:52:18.1 


-16:31 


:03 


10.46 


0.06 


4684 


2.92 


-0.21 


1.30 


1.26 


34 



Note. — Name, equatorial coordinates for epoch (J2000.0), B-magnitude mg, extinction A v , and derive 
atmospheric parameters are shown. Table 3 is published in its entirety in the electronic edition. 



Table 4. Mean stellar Parameters and uncertainties obtained from two-epoch observations. 



Name log g [Fe/H] £ vsini 



Tyc2773-01666-1 


2 


930 


± 





238 


-0 


279 


± 





021 


1 


583 


± 





091 


1 


68 


± 





16 


Tyc2269-01094-1 


3 


115 


± 





064 


-0 


405 


± 





031 





549 


± 





049 


3 


71 


± 





48 


Tyc2185-00133-1 


2 


774 


± 





340 


-0 


657 


± 





016 


1 


331 


± 





063 


5 


93 


± 





01 


Tycl726-00742-l 


2 


347 


± 





177 


-0 


294 


± 





042 


1 


333 


± 





001 


1 


94 


± 





20 


Tycl780-00654-l 


2 


444 


± 





354 


-0 


436 


± 





111 





953 


± 





248 


4 


63 


± 





07 


Tyc0001-00818-1 


2 


943 


± 





257 





161 


± 





060 


1 


228 


± 





210 





00 


± 





00 


G0410+11.148 


1 


785 


± 





327 


-0 


324 


± 





077 


1 


250 


± 





057 


3 


12 


± 





29 


G1705+22. 14898 


2 


267 


± 





299 


-1 


212 


± 





049 


1 


895 


± 





048 


5 


07 


± 





36 


G2238+33.1387 


2 


940 


± 





058 





132 


± 





044 


1 


031 


± 





088 





10 


± 





10 


Tycl753-00655-l 


1 


752 


± 





363 


-0 


319 


± 





067 


1 


385 


± 





034 


1 


50 


± 





18 


Tyc3973-02373-l 


2 


151 


± 





145 


-0 


232 


± 





056 


1 


411 


± 





093 


1 


39 


± 





04 


Mean 











238 











050 











089 











17 
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Table 5. Comparison of derived parameters for published values of iron abundance [Fe/H], 
surface gravity log g, and microturbulence £ to those obtained with TAGRA in this paper. 



Star 


T eff publ. 


[Fe/H] publ. 


log g publ. 


£ publ 


[Fe/H] 


log g 






K 


dex 


dex 


km/s 


dex 


dex 


km/s 


Sun 


5778 


0.00 


4.44 


1.10 


-0.06 


4.49 


1.09 


Arcturus 


4280 


-0.60 


1.90 


1.60 


-0.69 


1.78 


1.76 


7 Dra 


3830 


-0.14 


1.55 




-0.15 


1.25 


1.62 


11 Umi 


4120 


0.07 


2.03 




-0.23 


1.55 


1.69 


P Oph 


4475 


0.00 


1.70 




0.02 


2.33 


1.53 


HD152812 


4150 


-0.44 


2.10 




-0.62 


1.40 


1.56 


HD160290 


4440 


-0.21 


2.59 




-0.46 


1.84 


1.44 


HD218029 


4280 


0.07 


2.28 




-0.04 


2.04 


1.58 


7r Her 


4100 


-0.18 


1.68 




-0.11 


1.59 


1.84 


v Dra 


4520 


-0.12 


2.55 




-0.17 


2.18 


1.53 



Note. - Name of the star, its effetive temperature, and published [Fe/H], log g, 
and £ versus those obtained with TAGRA. 
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Finally, the accuracy of the TAGRA pipeline was tested from an analysis of the spec- 
trum of Arcturus degraded with the addition of artificial noise. We performed Monte Carlo 
simulations and derived the basic stellar parameters for spectra with a range in S/N from 10 
to 100, and produced sets of thirty noisy artificial spectra for various values of S/N. Figure @] 
shows the difference between the parameters log g, [Fe/H], £, and vsini estimated from 
the original high-S/N spectra and the same spectra with added noise. The solid line with 
squares shows the difference (zero at the vertical axes in all panels in Figure H] is designated 
by the dash-dotted line and corresponds to the case of no difference) . Dotted lines designate 
the r.m.s. scatter of the parameters determined from the degraded spectra. As can be seen 
from Figure HI the spectrum quality imparts rather small uncertainties in the metallicity, 
microturbulence and vsini (0.05 dex, 0.05 dex, and 0.2 km/s, respectively), with somewhat 
larger uncertainties in the surface gravity (0.25 dex). The vsini and log g also exhibit some 
moderate systematic differences for the degraded spectra (of the order of 0.1 km s _1 and 
+0.1 dex, respectively). 

Considering the reliability of basic atmospheric parameters derived using the TAGRA 
pipeline, we conclude that the metallicity ([Fe/H]), microturbulence (£), and vsini are the 
best-defined parameters. The surface gravities have poorer accuracies and may be affected 
somewhat by systematic errors. This is not unexpected because the primary source of log g 
is the balance between the Fel and Fell lines. There are not as many Fell lines as Fel in 
this wavelength range and, in addition, the available Fell lines are concentrated towards the 
more line-crowded blue end of the spectra. 



3.4. Magnesium Abundances 



The high-resolution spectra analyzed in this study can be used to probe chemical evo- 
lution using a number of different elements, and an analysis of Mg abundances, which is a 
product of Supenovae Type II (SN II), is presented here as a comparison to the Fe abun- 
dances, produced mainly in SN la. Within the spectral interval covered here, the Mg I line 
at 5711.090 A is well-defined and me asurable. Th i s line has a relatively accurate oscillator 
strength (~14%) and was used in the ISmith et al.l ( 120001 ) study of uj Cen. In the abundance 
analysis of Mg, the atmospheric parameters defined by the analysis of the Fe lines were used 
while the Mg abundance was varied to determine the best-fit profile to the Mg I line. 
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4. Discussion 

Figure [5] shows the distribution of T e ff, log g, and [Fe/H] for all sample stars. Since this 
study is comprised of two different samples, the GGSS candidates are indicated as the solid 
line histogram, while the Tycho stars are represented as the dashed line histogram. The top 
panel of Figure 5 shows clearly that the Tycho sample stars, which were picked to be probable 
clump red giants, tend to cluster in a narrower T e fj range, as would be expected for clump 
giants. The GGSS giants, on the other hand, spread across the entire red giant branch, as 
evidenced by the broad distribution in T e g. The middle panel shows the distribution of the 
samples in terms of their surface gravities. The average log g for the Tycho sample is 2.56 
± 0.55 and a similar average value is found for the GGSS sample 2.42 ± 0.64. The bottom 
panel illustrates the [Fe/H] distributions and shows the mean subsolar metallicities of both 
Tycho and GGSS samples, with these samples dominated by thick disk members. 

Figure [6] places the sample stars on a diagram of T e ff versus log g. Most of the stars are 
located along typi cal red giant b r anch and clump loci. The lower panel in F igure E] confirms 



one c onclusion by lBizyaev et al.l (120061 ) that the photometric GGSS survey (IPatterson et al. 



20011 ) selected red giants very efficiently. The Tycho sample was selected with color cuts and 



has no stars hotter than 4750 K in our sample of stars. The Tycho obj ects tend to concent rate 



in the location of the red giant clump in the T - log g diagram (see iBizyaev et al.ll2006l . for 
example). 

In Figure 7 we plot the values of [Mg/Fe] versus [Fe/H] for the high-S/N sample of the 
brighter giants and note that the well-established increase of [Mg/Fe] with decreasing [Fe/H] 
is recovered in this automated analysis. As an illustration of the agreement between this 
study and previous work, the abundances from Mishenina et al. (2006), who analyzed nearby 
bright clump giants, are also included in Figure 7. The general agreement demonstrates that 
the automated analysis techniques developed here produce reasonable abundance results, at 
least from spectra of relatively high S/N. 

Of more interest to continuing abundance studies based on these spectra are the Mg 
abundances derived from the entire sample, which generally consists of spectra with lower 
values of S/N. The abundance values of [Mg/Fe] versus [Fe/H] for all of the GGSS and Tycho 
giants are shown in Figure 8; here again, the well-established trend of increasing [Mg/Fe] 
with decreasing [Fe/H] is clearly evident. As a comparison study which is based on high- 
S/N spectra and employs a classical analysis, we include in Figure 8 the abundances from 
Reddy et al. (2003; 2006), which sampled thick and thin disk stars. This is an appropriate 
comparison since the majority of giants in this study are likely to be members of the thick 
disk. The summary point from Figures 7 and 8 is that this set of spectra can be used to 
probe the chemical evolution of a number of elements. 
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5. Conclusions 



We have developed a pipeline to estimate the stellar atmospheric parameters, metallic- 
ities and magnesium abundances for a large set of spectra of candidates to the SIM Planet 
Quest reference frame stars observed with an echelle spectrograph in 2004-2006. Most of 
the stars photometrically preselected as red giants turn out to be bona fide red giants with 
slightly subsolar metallicities, on average. The Tycho sample stars were observed with higher 
S/N, on average, and are more concentrated to the loci of red giant clump stars than the 
GGSS sample. The latter one consists of giant stars with T e ff uniformly distributed between 
4000 K to 5000 K. The iron abundances of the GGSS and Tycho samples are similar, with 
average values of [Fe/H] = -0.30 and -0.24, respectively. This average is more metal-poor 
than the value for the thin disk. 

The values of [Mg/Fe] versus [Fe/H] in the GGSS and Tycho program stars derived from 
the automated analysis d eveloped here when compared to results from the classical analysis in 



Reddy et al.l ( 120031 . 120061 ) is very good, demonstrating that the automated analysis produces 
reliable abundances for Fe and Mg. The next step in the automated analysis of this spectral 
dataset, which will be the topic of a future paper, will be abundance determinations of 
additional elements, with the particular elements chosen in order to study chemical evolution 
in this sample of predominantly thick and thin disk giants. 



We are grateful to our collaborators who helped with observations: William Cohran, 
Mike Endl, Nairn Baliber (all from UT Austin), and Nataly Petrova (McDonald Observa- 
tory). The GGSS follow-up observations presented here were supported by the JPL and 
NASA via grant 99-04-OSS-058. We are grateful to the staff of McDonald observatory for 
their help and support in photometric and high resolution spectroscopic observations. This 
publication makes use of data products from the Two Micron All Sky Survey, which is a joint 
project of the University of Massachusetts and the IPAC/Caltech, funded by the NASA and 
NSF. 
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Fig. I. — An echelle order of the Arcturus spectrum and the polynomial fitting defining the 
continuum (solid line). 
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Fig. 2. — Observed and best-fitting synthetic profiles for a sample of Fe lines in Tyc2630- 
01070-1, a typical star from our sample. The squares indicate the observed spectra, the solid 
curve shows the best-fit synthetic model, and the dashed curves are offset by ± 0.5 dex in 
[Fe/H] . 
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Fig. 3. — Comparison between the derived surface gravities (log g), metallicities ([Fe/H]) 
and microturbulent velocities (£) with published values from the literature for 10 bright 
stars which were observed as a control sample. The literature values are from Souburan et 
al. (2008), Luck & Challener (1995), McWilliam (1990) and Brown et al. (1989). 
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Fig. 4. — Differences between the parameters log g, [Fe/H], £, and vsini estimated from the 
noise-added spectra of Arcturus and those from high-S/N spectra (see explanations in text). 
The solid line with squares shows the difference. Zero at the vertical axes on all panels in 
Figure H] is designated by the dash-dotted horizontal line. Dotted lines designate the r.m.s. 
scatter of the parameters determined from the deteriorated spectra. 
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Fig. 5. — Distributions of T e ff, log g, and [Fe/H] for the GGSS (solid line) and Tycho 
(dashed line) samples. In terms of the average effective temperatures and surface gravities, 
the two samples are similar (< T e ff >=4567 and 4503 with < logg >=2.56 and 2.42, for 
Tycho and GGSS, respectively). The effective temperatures of the stars in the Tycho sample, 
however, are more localized as they are predominantly red clump stars. 
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Fig. 6. — The diagram T e ff versus log g for the Tycho (upper panel) and GGSS (lower 
panel) samples. The distribution of the two different samples in this diagram illustrates 
the clustering of the Tycho sample to a more localized region in the T e g-log g plane when 
compared to the GGSS sample. The Tycho sample contains a sizable fraction of red giant 
clump stars, while the GGSS sample consists of giants distributed more uniformly across the 
first ascent of the red giant branch. 
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Fig. 7. — The ratio of Magnesium/Iron abundances versus iron abundances in the control 
sample of 10 bright stars (filled circles), with the increase in [Mg/Fe] with decreasing [Fe/H] 
clearly evident. Included in this figure are abundance results from Mishenina et al. (2006) 
shown with diamonds and the agreement between the two studies is excellent. 
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Fig. 8.— The values of [Mg/Fe] versus [Fe/H] in the GGSS and Tycho program star s 
(dots). The squares and triangles show the abundance results from lReddy et al.l ( I2003l . l2006l ). 
respectively. The agreement between the abundances derived from the automated analy sis 
developed here when compared to the classical analysis from iReddy et al.l ( 120031 . 120061 ) is 
very good, demonstrating that the automated analysis produces reliable abundances for Fe 
and Mg. 



